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BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0003] The present invention relates generally to a system and method for 

performing noise cancellation in narrowband satellite communication systems. More 
particularly, the present invention relates to performing noise cancellation for estimated 
signal parameters on a demodulated signal. 

2. Description of the Related Art 

[0004] Efficient use of available bandwidth in wireless, such as 

satellite, communications applications is a problem of paramount 
importance. An example of such a narrow band satellite includes very 
small aperture terminal (VSAT) systems. VSAT systems use compact 
earth stations that are installed at one or more customer's premises to 
provide links among the premises over a wide coverage area. Typically, in such 



systems, remote ground terminals are used for communicating via a 
geosynchronous satellite from a remote location to a central hub station or 
other remote locations. The central hub station communicates with multiple 
remote ground terminals. VSAT systems are used to handle customer 
network requirements, from small retail sites up to major regional offices, 
and can support two-way data, voice, multi-media, and other types of data. 
A particular advantage of these systems is their relatively low site cost and 
small earth-station size. 

[0005] In wireless systems, multiple users share the same bandwidth. 

Charmel sharing through fixed-allocation, demand assigned or random- 
allocation modes is known as multiple access. Two of the more commonly 
known basic multiple-access techniques include time division multiple 
access (TDM A) and code division multiple access (CDMA). 
[0006] VSAT type systems have traditionally implemented TDMA 

using time division multiplexed (TDM) mode. Such systems generally are 
used for low speed (300 bps to 19,200 bps) data communications such as 
credit card processing and verification, point-of-sale inventory control and 
general business data connectivity. A typical TDM/ TDMA network, when 
implemented in a star topology (FIG. 1), uses a large satellite hub system 
that manages all network terminal access and routing. Data is transmitted 
to and from the hub in short bursts on satellite channels that are shared with 
a ntmiber of other VSAT terminals. The hub communicates with these 
VSAT terminals over a higher speed outbound TDM satellite carrier. The 
terminals transmit back to the hub on assigned inboimd carriers using TDM 
protocols. Such a combination enables a predetermined number of slots in 
time each second that each terminal can transmit data. In addition, more or 
less time can dynamicaEy be assigned to the terminals based upon each 
terminal's individual requirements. 

[0007] In contrast, in a CDMA type system a user's station signal is 
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multiplied by a unique spreading code at a high speed to be spread in a 
wide frequency band. Thereafter, the signal is transmitted to a transmission 
path. In a receiving side, the signal that was multiplexed by the spreading 
code is subjected to a despreading process to detect a desired signal. Signal 
detection is based on a unique spreading code assigned to a user's station. If 
despreading is carried out with reference to a particular code used to spread 
a transmission signal, a user's station signal is correctly reproduced. 
[0008] Regardless of the access technique used, increased efficiency 

and lower cost is a primary goal. Accordingly, efficiencies in bandwidth 
may be realized using techniques such as crowding of adjacent channels, 
frequency re-use, and increasing of data rates, generally resulting in an 
increased amount of data traveling through the limited amount of available 
bandwidth. Unfortunately, however, such techniques introduce a 
sigiuficant amotmt of interference which must be canceled. Combined 
multi-user detection and decoding is believed to have the potential to 
improve performance to match that of an interference-free system. 
However, most development heretofore has been in CDMA based systems. 
In particular, it is known that the optimum receiver for CDMA system 
employing Forward Error Control (FEC) coding combines the trellises of 
both the multi-user detector and the FEC code. However, the complexity of 
such a receiver is exponential in the product of the ntimber of users and the 
constraint length of the code. This complexity makes the use of the optimal 
detector prohibitive for even small systems. 

[0009] Copending U.S. Patent Application entitled "Combined 

Interference Cancellation with FEC Decoding for High Spectral Efficiency 
Satellite Communicatioris" Serial No. 09/436,670, filed November 10, 1999 
discloses a method and apparatus for a low complexity cancellation scheme in 
narrow band type satellite applicatioiis that allows for efficient utilization of 



available bandwidth by eliminating the interference resulting from the 
aggressive channel crowding. 

[0010] However, there is a need for performing iterative noise 

cancellation on signal parameters such as frequency, timing phase and signal 
strength. The invention should preferably use a demodulator to achieve this 
result. 

SUMMARY OF THE INVENTION 

[0011] Briefly, the present invention relates to a satellite 

commtmications system and method for achieving efficient utilization of 
available bandwidth for satellite applications such as fixed wireless, mobile 
satellite systems and other narrow-band type applications. A soft decision- 
feedback scheme is used iteratively in combination with Forward Error 
Correction (EEC) decoding for interference cancellation to enable efficient 
use of the available bandwidth through aggressive crowding of adjacent 
charmels. 

[0012] In a first embodiment of the present invention, a multiple 

channel decoding receiver is provided. The multiple channel receiver 
includes a matched-filter bank that is used to receive signals and provide 
initial estimates of data. The estimates are fed to an interference canceler, 
which provides channel information to a demodulator. The demodulator 
outputs a signal representative of an estimation of at least one parameter of 
the channel information. The estimated signal is then provided to a decoder 
which calculates an estimated interference value based on the estimation 
signal. The estimated interference value comprises soft interference 
estimates which are fed back into the interference canceler. The soft 
estimates of the interfering signals are subtracted from the matched-filter 
outputs in the interference canceler to generate new, refined, soft approxi- 
mations of the data. These refined soft estimates are then fed to the chaimel 
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demodulators for the next iteration. The process is repeated iteratively until 
the channel interference reaches acceptable levels. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is pointed out with particularity in the 
appended claims. However, other features of the invention will become 
more apparent, and the invention will be best understood by referring to the 
following detailed description in conjunction with the accompanying 
drawings, in which: 

FIG. 1 illustrates a VSAT system in a star topology; 

FIG. 2 is a schematic diagram of an exemplary ACI transmitter 
model in accordance with an embodiment of the present invention; 

FIG. 3 is a block diagram of a first exemplary embodiment of 
an ACI receiver model in accordance with the present invention; 

FIG. 4 illustrates a spectral view of a signal model for known 

systems; 

FIG. 5 illustrates a spectral view of a signal model in 
accordance with an embodiment of the present invention; 

FIG. 6 is a block diagram of a multi-channel receiver which 
combines interference cancellation with forward error correction (FEC) 
decoding in accordance with the present invention; 

FIG. 7 is a block diagram of a single-channel receiver which 
combines interference cancellation with forward error correction (FEC) 
decoding in accordance with the present invention; 

FIGS. 8a and 8b are graphs illustrating the performance of the 
proposed iterative multi-charmel receiver at channel spacing of 0.757;"^ with 
4 states and 16 states convolutional code, respectively; and 

FIG. 9 is a graph illustrating the performance of the proposed 
iterative single-charmel receiver at charmel spacing of 0.757;""^ with 4 states 



convolutional code. 



DETAILED DESCRIPTION OF THE INVENTION 
[0013] The present invention relates to a satellite commxinications 

system and method for achieving efficient utilization of available bandwidth 
for satellite applications. In particular, a soft decision-feedback scheme is 
used iteratively for interference cancellation in combination with FEC 
decoding to enable efficient use of the available bandwidth using techniques 
such as crowding of adjacent channels, frequency re-use, and increasing the 
data rates. A particular advantage of such a system is the ability to 
eliminate interference, such as adjacent channel interference (ACI), that may 
be introduced during, for example, channel crowding, thereby resulting in a 
higher spectral efficiency. For example, the present embodiment enables a 
satellite system to operate at a bandwidth efficiency level of 2.66 bits-per- 
second/Hz with minimum additional energy requirement in the signal-to- 
noise ratio range of interest using only a four state FEC code. This 
corresponds to an approximately 55% improvement in spectral utilization 
over current systems that employ similar modulation techniques. This 
improvement is expected to increase when more efficient FEC codes are 
used. 

[0014] Turning now to FIG. 2, there is illustrated an exemplary VSAT 

system communicating in a TDMA mode. It is to be noted that the 
principals embodied in the present embodiment may also be applied to 
other satellite communications systems as well. As shown in the figure, the 
VSAT system, such as available from Hughes Network Systems, includes a 
central hub station 102 that controls one or more earth stations 104A-104B 
located on customers' premises. The earth stations 104A-104B and the 
central hub station 102 communicate with each other using a 
geosynchronous satellite 106. Each of tiie earth stations 104A-104B has a 



receiver 108A-108B for receiving and decoding signals received from the 
satellite 106 and transmitters llOA-llOB for transmitting data to the satellite 
106. The hub, or base station, station 102 similarly includes a receiver 112 
for receiving and decoding signals received from the satellite 106 and a 
transmitter 114 for transmitting data to the satellite 106. 
[0015] Turning to FIG. 3, an exemplary ACT transmitter model 120 is 

shown which may be used in the earth station transmitters llOA-llOB and 
the hub transmitter 114. The transmitter 120 receives data from a first 
source 122A to an Mth source 122C. Converters 124A-124C convert the data 
from a binary phase shift keying (BPSK) signal to a quadrature phase shift 
keying (QPSK) signal. The resultant frequency domain pulse 126 is 
interleaved and transmitted as signal s(t) 128. The signal s(t) 128 which 
models the situation of adjacent chaimel interference caused by signal 
crowding, consists of the signal in noise as 

m = s(t) + n(t) (1) 

[0016] The n(t) is the standard additive white Gaussian noise 

(AWGN) with single-sided power spectral density level of N^, (Watts/Hz). 
The signals s(t) models the situation of ACI in which there are M adjacent 
data sources that are identical and independent. Each source transmits a 
QPSK signal at the rate of T;^ vdth an arbitrary unit-energy pulse shape p(t). 
The signal is described in complex form as 

s(t)^ Re{s(Oe'^"-^ci} (2) 
where/^ is the carrier frequency and s(t) is the baseband complex envelope 
of the signal and is mathematically expressed as 
M 00 

m=l A:=-oo 

The data streams {a^j; m=l,...,M} consist of QPSK symbols taking on the 
equi-likely values (±1 ±j) and are statistically related as 



^i«*«.,/««jj = |0, otherwise (4) 

[0017] Assuming without loss of generality that M is an odd integer 

and that the desired center, or (M + l)/2-th, conveys the desired data and 
that the other signals are viewed as being adjacent-channel interferers ((M- 
l)/2 ones on either side), namely A = <...< Af (^.i)/2<...< Af^,. The present 
embodiment is used when the channel spacing values, compared with that 
of the center channel, are small enough to cause large amounts of overlap in 
the spectra. Furthermore, in practical systems, these channels are equally 
spaced in frequency, say by Af . In terms of the above, then 
, ^ ( 71/ + ll 

^fm = {^-—Y~)^f'^ m=l,2..;M (5) 

[0018] It is to be noted that M should not be interpreted as being the 

number of channels in the entire available bandwidth. Instead, it is the 
number of channels that the receiver wishes to process jointly to armounce a 
decision regarding the desked data stream. Although the use of guard 
bands is known when separating channels, their use consumes a non-trivial 
amount of bandwidth, thereby decreasing spectral efficiency. Therefore, in 
the absence of an installed guard band, the outermost or "edge" channels, 
i.e., 1st and Mth, will always have interference. However, as described 
further below, the present embodiment does not require that the receiver 
compensate for these edge charmels. For example, M is chosen to be seven. 
As such, the receiver processes seven charmels in the presence of two 
additional signals. Note that in the presence of guard band at the edges of 
the M charmels, the receiver will jointly receive all the M channels. This 
feature is particularly useful for the base-station which is interested in 
decoding more than one channel. 

[0019] Two criteria of performance are considered. The first is the 
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bandwidth efficiency, r|, in bits-per-second/Hz defined as the ratio of the 
data rate to the bandwidth used. This quantity is obtained in this case as 



[0020] The second is the energy efficiency defined as the signal-to- 

noise ration per bit required to achieve a specific bit error probability Pb(E) 
of the desired channel. Other measures of performance may also be used, 
such as symbol error probability and word error probability. 
[0021] It can be seen that the interference-free performance is 



(7) 



where E,, is the average bit energy and Q(x) is defined as the Gaussian 
probability integral 



<» 1 

'r^ X 

The performance described by the above equations is used as a benchmark 
to quantify the ability of the proposed receivers to suppress interference. 
[0022] From statistical theory of hypothesis testing, it is known that 

the optimum receiver is the one that minimizes sequence error probability 
and is derived from implementing the average likelihood-ration function 
(ALF). The likelihood function conditional on knowing a given signal is 
AWGN is then obtained as 



.(r(0)=exp|--^ J \rit)-7{tfdt 



where r(t) is the baseband complex envelope of the received waveform. 



(9) 
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Thus, 



ln(A(r(/))) 



(10) 



±1 



^(0 



M 

m=\ k-- 



i= \ k=-m 



By absorbing terms that are independent of the sought symbols and after 
some algebraic manipulations, maximizing the ALF over the information 
symbols is equivalent to maximizing the quantity J(ai,a2,...aJ of the metric, 
defined as 



MRe 



k=-co 



- z z 

A;=-oo /=-« 



M 

m=l 
M M 



(11) 



The above expression shows that {x^((k+e)T,); m=l,...,M} is a set of sufficient 
statistics which consists of an exemplary bank of matched filters 132, as 
shown in FIGs. 4 and 5, matched to the modulating signal in each channel, 
then sampled at the symbol rate of T^'^ . More generally defined as 
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(12) 



and also 



(13) 



The function of the optimal rule, or the maximum likelihood sequence 
estimation receiver 136, is to determine the sequence of information symbols 
(ai,a2,...aj^) that maximizes the metric shown above. If there are N symbols 
in a frame, then the most straightforward way of implementing the 
optimum receiver requires 4^ computations of the metric. However, this 
procedure can be implemented in the most efficient way by generalizing the 
modified Viterbi Algorithm (VA) of G. Ungerboeck, "Adaptive Maximum- 
Likelihood Receiver for Carrier-Modulated Data Transmission Systems," 
IEEE Transactions on Communications, pp. 624-636, May 1974. 
Accordingly, the metric shown above can be made recursive by the relation 

(14) 

Jp{a\^p,a2^p,...,aMp) 
= Jp- 1 {ai^p^ 1 , a2^p- 1 , . . . ,aj^ p_ j ) 

M 

^^m(iP^^m)Ts)- S Cn,mi'<P^ ^n)TsUp+ s^)Ts)an^p 



+ Re{ Z 



M 



1 ^n,m(ik^ Sri)T,^{p^ s^)T,)a„^j, 1 



Equivalentiy, the interference channel whose impulse response spans L 
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symbols can be viewed as a finite-state discrete-time machine where the 
state at discrete time i is defined as 

The VA then tracks the paths through the treUis and provides the solution to 
the problem of maximum-likelihood estimate of the state sequence. Thus, it 
is clear that the trellis has a maximirai of 4'^ states. Note that the efficiency 
of this modified VA stems from the fact that maximizing the likelihood 
function requires computing N4^ instead of 4^^ metrics, wherein L is 
typically much smaller than N. Thus, even-though the complexity remains 
exponential in the number of channels, making the optimum rule 
computationally intensive, the complexity becomes independent of N. It is 
to be noted that reduced-complexity versions of the vector VA, which use 
decision-feedback on a per-survivor basis, may also be used. It is worth 
noting that this receiver does not include the FEC decoding which is done 
separately. This will result in a loss in performance. The reason for this is 
the huge complexity of the optimum joint receiver. In this invention we 
present a low complexity receiver for joint demodulation and charmel 
decoding that achieves very close performance to the optimum receiver. 
[0023] Before providing a suboptimal cancellation structures, the 

effective channel that appears at the output of the matched-filter bank, 
{Xj„(t);m=l,2,...,M}, is characterized to yield 

M 00 
.n=l /=-oo 



From the above equation it is clear that the equivalent lowpass interference 
channel is described by the previously defined impulse response Cy(ti,t2), 
which can be equivalently represented in the frequency domain as 
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Q,/(^/-^2)a[ j P(/)P* (/- (A/,- - A/^))e^2-/^^2-0V/j 
xexp(-y(2KA/y - A/;)f2 + (^y - ^/)) 

This function represents the effective channel impulse response at the output 
of the jth matched-filter when excited by the Ith data source. It consists of 
the cascade of the pulse-shaping filter and the complex multiplier at the 
transmitter side, the channel, and the matched-filter at the receiver. It is to 
be noted that the impulse response in this case is time-varying, a condition 
that results from the presence of complex-exponential multipliers (or 
frequency shifters) in the system. As the channel spacing is increased, the 
magnitude of the impulse response decreases but its duration is increased, 
resulting in an equivalent channel v>rith larger memory span. 
From above, 

nm(0^jn(a)p* ia- t)e~-''^^"^^'""^^'"^da (18) 

wherein, the {n„,(t);m=l,2,...,M} is a set of zero-mean complex Gaussian 
random processes with covariance 

E{n'^l(t)nj(f)} = NoCij(t,f) (19) 

Several asstimptions are made to simplify notation. First, in regards to the 
piilse shaping, it is assimied that the pulse p(t) selected satisfies the Nyquist 
criterion of zero inter-symbol interference. This criterion is expressed in 
time as 

00 f 1 0 

-00 

or in frequency as 
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(21) 



This states that the aliased or folded version of the auto-power spectrum 
associated with the selected pulse must be flat. (The aliased version is what 
results when replicating the function at multiples of the symbol rate T/\) 
Note that the root-raised cosine pulse, which is a practical and bandwidth- 
efficient choice, is defined in the frequency domain as 

1- 0<|/|< (1-^)727; 

(1 - ^) / < (1 + 
0, otherwise 



Pif)- 



)/27; 



(22) 



where p is the roll-off parameter. Second, the spectral overlap of these 
channels does not exceed 50%. This, along with practical values of the roll- 
off parameter, indicates that the ACI on a given channel results from one 
adjacent interferer on either side. Third, it is assumed that the ACI extends 
over a finite time interval spaiming L symbols. The actual value of L is 
directly related to the amount of spectral overlap that exists between the 
channels. From basic principles of Fourier transforms, the value of L, which 
can be thought of as the memory of the interference channel, is larger for 
smaller overlap. 

[0024] Fourth, it is assumed that the receiver is able to maintain phase 

coherence and time synchronism. For the time synchronism situation, a 
situation representative of a satellite down-link application, it is assumed 
that the relative time delays are zero. Based on the above set of 
assumptions, the matched-filter statistic is described as 
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L 

+ Y.^m-l,k-iCm-\,m{i.k-i)Ts>kTs) 

'^l^ (23) 
+ E«m+U-zQ«+l,m((^-07;'A:rj 



The first term on the right-hand side of the above equation is the desired 
information symbol; the second term is the ACI contribution from the left 
charmel; while the third term is the ACI contribution from the right channel. 
The ACI is determined by the symbol-spaced samples of the cross- 
correlation between transmit and receive filters. The set 
{n^(kTs);m=l,2,...,M} has elements that are statistically correlated across 
different m's but independent for a specific m. 

[0025] Combined multi-user detection and decoding has received 

considerable attention recently with its potential to improve performance to 
match that of the single-user system. However, most of the work done has 
focused on spread-spectrum CDMA systems. It has been shown that the 
optimum receiver, for a CDMA system employing FEC coding, combines 
the trellises of both the multi-user detector and the FEC code. The 
complexity of this receiver is exponential in the product of the number of 
users and the constraint length of the code, making the implementation of 
the optimal detector prohibitive for even small systems. 
[0026] FIG. 4 which illustrates a spectral view of a signal model for 

known systems highlights the inefficiencies of the ctirrent art. As previously 
discussed above, the bandwidth is limited and an ideal scenario is to pack as 
many channels as possible into the limited available bandwidth. However, 
limitations such as interference between channels limit the number of 
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channels that can be packed into the available bandwidth. 
[0027] FIG. 5 illustrates a spectral view of a signal model in 

accordance with an embodiment of the present invention. Specifically, FIG. 5 
depicts a 40% improvement in bandwidth efficiency over prior art systems. 
Channels are packed closer together in the limited available bandwidth 
resulting in greater bandwidth capacity. By decoding additional channels, 
an even greater improvement can be achieved. 

[0028] Recently, a new powerful class of concatenated convolutional 

codes was proposed which use parallel concatenation of two (or more) 
recursive systematic convolutional codes (constituent codes) fed by two 
information sequences, of which the second is obtained from the first 
through the interposition of a long interleaver. One of the key factors 
contributing to the remarkable performance of this coding scheme is the 
elegant iterative Soft-Input Soft-Output (SISO) decoding stiructure whose 
performance was shown, via simulation, to approach that of the maximum- 
likelihood decoding, at signal-to-noise ratios very close to the Shannon limit, 
with much less complexity. This decoder is based on iteratively decoding 
the component codes and passing the so-called exti-insic information, which 
is a part of the component decoder soft output, to the next decoding stage. 
The impressive performance, achieved by this iterative decoding 
architecture, has encouraged several researchers to consider applying the 
same principle in the other sub-modules of the receivers. 
[0029] In this section, several strategies for joint iterative interference 

cancellation and FEC decoding, which provide a tradeoff between 
performance and complexity, are investigated. We will restiict otirselves to 
binary convolutional codes; the extension to higher order codes is 
stiraightforward. The proposed algorithms function independentiy of the 
type of channel decoders in the sense tiiat any SISO decoders can be used 
(for example, MAP, Log-MAP, or SOVA). The choice of a particular decoder 
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is determined by the allowable complexity at the receiver. Based on the 
Gray mapping assumption, each QPSK symbol is 

^m,k = ai^k + ja^,k (24) 

where aj„ ^. , and a^ j^ e {- l,l} • Without loss of generality, we will derive 
the different detection rules for the in-phase A:-th binary symbol associated 
with the m-th channel. 

Denote, 



The set A' can be similarly defined. Using the independence assumption, 
which is justified by the interleaving, the maximum-a-posteriori (MAP) 
detector is given by 

(26) 





^(^, 


n,k 




v-\,m+l},k-L<l<k+L A'^j,l)P\ 








n,k 




n-l,m+\],k-L<l<k+L A^j,l)^[ 





where is the updated log-likelihood ratio, P(Zm,k\^ is the 
conditional Gaussian distribution of the matched filter output as per (26). In 
(26), P((Zy /) and P(a^^) are obtained from the soft outputs of the 
previous iteration as follows 

^(^y,/ = ^J,l = na],i)P{aQj) (27) 



^(4/=!)=-^ (28) 
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l+e j'^ 



(29) 



where Aj i is the output log-likelihood ratio of the previous iteration. It is 
clear that results similar to (28), (29) hold for P{a9^) . 
[0030] The MAP detector requires a complexity of the order 

0(2 X 4-^(^^'^^^ ) which can be prohibitive for practical applications. There- 
fore, in the following, a lower complexity detection rule based on the MMSE 
principle is developed. The iterative MAP detection rule has been proposed 
for CDMA signals. 

[0031] This scheme, depicted in FIG. 6, uses the soft information 

supplied by the M single-user decoders to calculate the optimum, feed- 
forward and feed-back, filter weights after each iteration. FIG. 6 shows the 
multi-channel receiver 130 which combines interference cancellation with 
forward error cancellation (FEC) decoding using a matched filter 132 and 
channel estimation 136 as inputs to an interference canceler (IC) 134. 
Multiple channels M are provided from the interference canceler 134 to 
demodulators 138a... 138b which provide outputs to optional equalizers 
140a.. .140b. The equalizers 140a.. .140b provide outputs to deinterleaver 
142a... 142b, which provide outputs to SISO decoders 144a... 144b and 
interleavers 146a...l46b for outputs 1 through M. The interleaver output 
146a and 146b are fed back to the interference canceler 134 in a closed loop. 
[0032] The received signal is processed at the match filter 132 where 

interference associated with the received signal is suppressed or the 
transmitted pulse of the signal is matched. Specifically, the out of band 
interference is suppressed. The suppressed signal is then provided to the IC 
134 where a portion of the interference is canceled. 

[0033] The suppressed signal is provided to the demodulators 138a 

through 138b, where timing, phase, frequency and signal strength 
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estimation is performed on the signal. The estimated parameters are used by 
the IC 134 to regenerate the interference signal. The interference suppressed 
signal is used at the successive iteration by the demodulators 138a, 138b to 
re-estimate the parameters again. With each iteration, the demodulators 

138a, 138b encounter less interference. 

[0034] The interference can not practically be eliminated completely. 

Each iteration reduces the amount of interference. However, there are trade 
offs. For each iteration, the reduction of the amoxmt of interference becomes 
less. In addition, different parameters can have different convergence rates. 
For example, the amplitude convergence can gain satisfactory accuracy in 
the first few iterations while the phase and timing estimation can require 
further iterations. To accommodate this situation the demodulators 138a, 
138b functions can be terminated earlier than that of the decoders 144a, 144b 
and IC's 134 functions. In other words, the estimation is discontinued at the 
demodulators 138a, 138b for the parameters that converge. 
[0035] In order to limit the number of iterations performed on 

adjacent channels, in one embodiment of the present invention, the multiple 
channels M are grouped into carrier groups. Within each carrier group, 
interference cancellation techniques are performed to reduce the spacing 
between channels within the carrier groups. However, sufficient spacing is 
utilized between carrier groups so that interference cancellation is not 
required. In other words, sufficient spacing is provided between the carrier 
groups so that minimal interference occurs between carrier groups. 
[0036] If the multiple channels M were not grouped into carrier 

groups, the possibility of continuous interference cancellation could occur. 
That is, in large systems having a number of channels, multiple occurrences 
of interference could occur. The interference cancellation process would 
continue the process for each channel. By having the multiple channels M 
grouped into carrier groups, the number of iterations is limited by the 
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carrier group size. 

[0037] In another embodiment of the present invention, each carrier 

group is further divided into two subgroups - an even and odd group. 
Demodulation/ decoding can be performed first in one subgroup. When the 
receiver processes the second sub group, the interference information gained 
from processing the first subgroup can be utilized for caricelling out part of 
the interfereiice in the second subgroup. In this manner, the second 
subgroup can benefit from an extra iteration without having to perform that 
extra iteration. 

[0038] Equalizers 140a, 140b are optional and equalize the received 

signal from demodulators 138a, 138b. In one embodiment of the invention, 
equalizers 140a, 140b can be used where there is severe inter-symbol channel 
interference. The equalizers 140a, 140b are depicted as being connected 
between demodulators 138a, 138b and decoders 144a, 144b. However, the 
equalizers 144a, 144b can be inserted after and/ or in parallel with the 
demodulators 138a, 138b and can be included with the IC 134. Thus, the 
equalizers 140a, 140b can benefit from the reduced interference level and the 
reliability information from the soft output decoder 144a, 144b. 
[0039] However, it is not just the equalizers 140a, 140b which benefit 

from the reduced interference level and the reliability information from the 
soft output of the decoders 144a, 144b, the demodulators 138a, 138b also 
benefit. The decoders 144a, 144b provide reliable hard decisions for the 
tracking of loops for the demodulators 138, 138b. For example, the decoders 
144a, 144b provide information on what's being transmitted e.g., a zero or 
one. Further, the decoders 144a, 144b also provides soft information on the 
probability of the information being correct e.g., the information being 
transmitted has a probability of .7 of being a one. This provides reliability to 
demodulators 138a, 138b. 

[0040] Although the present invention is discussed in terms of linear 
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systems, it will be appreciated by those skilled in the art that the present 
invention is applicable to nonlinear systems. The present invention can be 
applied to any FEC based transmitter system where the transmitter 
characteristics are known or can be estimated by the receiver. 
[0041] FIG. 7 on the other hand shows the single-channel 

embodiment of a receiver 148 which combines interference cancellation and 
FEC decoding employing the match filter bank 132 with a minimum means 
square error (MMSE) transversal filter 150 which provides a single charmel 
to a deinterleaver 154, SISO decoder 156, and interleaver 158 which provides 
input to calculate feed forward and feed bade coefficients at 160. The feed 
forward and feed back coefficients calculation 160 also receives the channel 
estimation 152 to provide feed forward and feed back coefficients to the 
MMSE transversal filter 150. For this we extend the algorithms for CDMA 
signals to the current narrow band TDMA application. Let J be a 
[m{2L + l) X l] complex vector of the matched filter bank outputs from the 
(k-L)-th to the (k+L)-ih samples. Then / can be written as 
M k+2L 

1=11 2 . lajj + n (30) 
j=l l=k-2L 

where, for example, 

Cj,M((^+ ej)T,,(k+L+ emWs) f 

is the ACI associated with a j i , and n is the Gaussian noise vector. Let, R^, 
the interference correlation matrix, be defined as 

and 
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(32) 



then 



Now, the updated decoder input is calculated from 
ym,k = £fZ + Cb 



(33) 

(34) 
(35) 



. = Im 



ym,k\ 



(36) 



where £y is the [ M{2L + 1) x l] feed-forward coefficients vector, cjy is the 
feed-back coefficient. Restricting the filter to have a single feed-back 
coefficient, rather than a vector, should not result in a loss of degrees of 

freedom. The coefficients £f , cjy are obtained through minimizing the 
MSE between the data symbol and its estimate, given by 

\ym,k-am,kf^ 

~\ rr. |2" 



MSE= E 
E 
E 



\£fl+%-am,k\ 



(37) 



It may be shown that the optimum solution based on MMSE must satisfy the 
following relations 



(38)- 
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^[{£}{c^,lc^,„,k + «)+ c, - a„^,)z_''] = 0 (39) 

Note that the relation (38) ensures that the output of the MMSE filter is 
unbiased, while (39) is a direct application of the orthogonality principle. 
Solving (38), (39), we obtain the following results for the feed-forward and 
feed-back filter coefficients 

= -c^R^Elq] (40) 

c} = C^j^[a CC^y ^ (41) 

where, by definition, we have 

^iQm,k^m,k (42) 

BAR^ E[ga^]R^^ (43) 

CAR^E[q] (44) 

and j^n is the [ M(2L + 1) x M(2L + 1)] noise covariance matrix which 
may be constructed using a component-wise relation. In (43)-(44), the Elg} 
and E[aa^] values are obtained from the following component-wise 
relations 

Eiajj)= E{ajj)^jEia%) (45) 

^(^il)=—l (46) 

e'-J'' + 1 

Eiajjajj) = 1 (47) 
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(48) 



where (48) follows from the independence assumption. In the first decoding 
iteration, we select £[0] = 0 . The feed-forward filter coefficients vector, 
£f , in this iteration is given by the MMSE equations derived in (41), and 
the feedback coefficient c^, = 0. After each iteration, the E\a\ values are 
re-calculated using the decoder's soft outputs. The E\q\ values are then 
used to generate the new set of filter coefficients as described. In the 
asymptotic case where \E\a\ = 1 , the receiver is equivalent to the 
subtractive interference canceler. This is expected, since \E\qX = 1 means 
that previous iteration decisions, for the interference, are error-free. In this 
case, the subtractive interference canceler is the optimum solution. The two 
main sources of complexity in this algorithm are the matrix inversion 
operation required in (41), and the need for M SISO channel decoders. In 
the following, we will investigate how to lower the computational 
complexity by proposing a soft interference cancellation algorithm that does 
not require a matrix inversion operation. In addition, the performance of 
the MMSE SISO detector will be studied assuming single-channel decoding. 
Based on (40), (41), yfn,k can be written as 

ym,k = £'^(Z-R^E[q]) (49) 

By observing that the matrix inversion operation is only required to 
compute c J , the following approximation is proposed 
T _ T 

Cf - j)(2z+ 1)+ z+ 1 (50) 
where e(^_iX2Z+l)+Z+l = [0.-»0,l(^_ix2i+i)+z,+i 0,...,0] . Hence, 

ym,k = Xm,k - iR^E{q\) (51) 

The complexity of this algorithm is a linear function of the product of the 
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number of interfering users and the interference memory (i.e., 0(2( L + 1)) ). 
This algorithm can be regarded as a soft subtractive interference canceler. 
This is so as the decoder's soft outputs are used to calculate estimates of the 
transmitted symbols, E[a] ; the estimates of the transmitted symbols and 
the interference cross-correlation matrix, are used to generate updated 
estimates of the interference signals, at the output of the matched filter. The 
interference estimates are then subtracted from the matched filter output 
^m,k resulting in the next decoding iteration input. It is interesting to note 
that this MMSE-based development results in a scheme similar to the one 
proposed under a different derivation for CDMA signals. 
[0042] Although the iterative SISO MMSE algorithm was developed 

assuming the use of M charmel decoders, it is straightforward to modify the 
algorithm to be used with M matched filters but with an arbitrary number of 
decoders that is smaller than M. For any undecoded user d, we have 

^l^dj 1 = 0 for any I (52) 

after any iteration. It is particularly interesting to consider the case of single- 
channel decoding which is displayed in FIG. 7. In this case, the soft 
information is fed-back to cancel the ISI, appearing at the matched filter 
bank outputs, and the feed-forward MMSE filter coefficients are used to 
suppress the ACL It is easy to see that, assuming error-free feedback, this 
algorithm is capable of suppressing the ACI, except for the edge chaimels 
effect, asymptotically. 

[0043] Monte-Carlo simulations are implemented to evaluate the bit 

error rate performance and demonstrate the effectiveness of the proposed 
solutions that combine interference cancellation and FEC decoding. The 
FEC encoding considered is the optimum rate V2 convolutional code with 4 
and 16 states and the decoder uses the soft output Viterbi algorithm (SOVA). 
The channel spacing used is A/ = 0.7527 ^ which corresponds to a 
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spectral efficiency level of 2.6671 bps/Hz or an improvement of 55% 
compared to the current state-of-the-art. 

[0044] FIGS. 8a and 8b compare the performance of the soft 

interference cancellation scheme, the conventional receiver, and the 
interference-free system. The receiver processes seven channels jointly in a 
presence of a total of nine QPSK sources. It is clear that the performance of 
the proposed iterative decoding and interference cancellation algorithm is 
better than the conventional receiver and very close to the interference-free 
system, with a difference of less than 0.5 dB using four iterations when the 
input SNR is about 4 dB. It is also noted that as the SNR increases, the 
performance gap between the iterative algorithm and the conventional 
receiver increases while it diminishes more compared with the interference- 
free system. 

[0045] FIG. 9 includes the performance achieved by the iterative 

MMSE algorithm, assuming single-channel coding, the conventional 
receiver, the feed-forward MMSE receiver, and the interference-free system. 
The receiver has three matched filters and it was assumed that only three 
channels are transmitting simultaneously (i.e., neglecting the edge effect). It 
is clear that the proposed algorithm provides considerable gain in 
performance compared to the conventional receiver and the feed-forward 
MMSE receiver. However, the difference in performance between the 
interference-free bound and the single-channel decoding algorithm is 
between 1.5-2 dBs which may be unacceptable in some cases. The 
performance can be improved upon by decoding more charmels, offering a 
tradeoff between performance quality and receiver computational load. 
[0046] The present invention and its performance have been 

described primarily in association with the customary AWGN channel. It is 
however clear from concepts in this patent to extend it to the case of fading 
channel. This can be done, for example, by estimating the chaimel directly 
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and incorporating this estimate into the reconstruction of the interference for 
subsequent cancellation. 

[0047] In order to prevent the synchronizing sequence from being 

obscured by the adjacent charmels, a system level mechanism is needed to 
ensure that only one user is transmitting during its known symbol period. 
[0048] Obviously, many modifications and variations of the present 

invention are possible in light of the above teachings. Thus, it is to be 
understood that, within the scope of the appended claims, the invention 
may be practiced otherwise than as specifically described above. 



